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I summarize recent lattice results on QCD at finite temperatures and densities. Studies 

on the nature of the QCD transition at the physical point, continuum extrapolations of 

thermodynamic quantities, and new calculations of hadronic spectral functions applying 

the maximum entropy method are discussed. Recent advances in finite density QCD are 

^ also reviewed. 
O 
(N 

O : 1. INTRODUCTION 

Lattice QCD provides us with the most systematic way to calculate non-perturbative 
hadronic quantities directly from the first principles of QCD, without resorting to models 
and effective theories. Because a simulation of the fully realistic case is still difficult, we 
have to carry out several extrapolations to extract physical quantities of phenomenological 
relevance. These include the continuum extrapolation, the extrapolation to the physical 
quark mass point, and, for heavy-ion collisions, the extrapolation to finite quark densities 
also. Recent progress in the computer power and simulation techniques have enabled us 
to put well-controlled systematic errors to several quantities. 

In this article, I review the status of lattice QCD at finite temperatures and densities, 
concentrating on studies which appeared after the Quark Matter 2001 conference [ |TJ. 
In Sec. 0, I summarize the developments in finite temperature QCD, selecting the topics 
of the approaches toward the physical point and the continuum limit. I also discuss 
recent studies of hadronic spectral functions adopting a new calculation technique, the 
maximal entropy method. A noticeable trend of lattice studies in the last 1.5 years is a 
rapid increase of investigations at finite densities, induced by a study of QCD at small 
chemical potentials by Fodor and Katz [ |2|. I introduce these developments in Sec. |3|. 
Short conclusions are given in Sec. B. 



2. FINITE TEMPERATURE QCD 

2.1. Toward the physical point 

Because the cost to simulate dynamical quarks is quite high on the lattice, most large 
scale simulations are done in the approximation that only degenerate u and d quarks 
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Figure 1. Nature of the QCD transition as a function of quark masses, (a) Theoretical 
expectations from various effective models, (b) Phase structure around the three flavor 
critical point from a staggered quark simulation using the standard action on an N t = 4 
lattice f [T3 



are treated dynamically (Np = 2 QCD). However, since m s ~ 100 MeV is around the 
transition temperature T c ~ 150 MeV, inclusion of the dynamical s quark is important 
(Np = 2 + 1 QCD). Fig. [L](a) summarizes our expectations on the nature of the QCD 
transition as a function of quark masses based on various theoretical studies of effective 
models. 

In the limit of heavy s quark (Np = 2 QCD), the chiral transition at m uc i = is expected 
to be second order in the universality class of the three dimensional 0(4) Heisenberg 
model [ [|. Previous studies of two- flavor QCD have been done mostly on N t = 4 and 
6 lattices. Through improvement of the lattice theory, several quantities, such as the 
transition temperature and the EOS, have started to show convergence among different 
lattice formulations [ f|, On the other hand, the expected 0(4) scaling could be seen 
only with improved Wilson- type quarks [ @, U, and so far not with staggered quarks [ 

When we decrease m s along the m u d = line, we will encounter a tri-critical point m* 
where the second order chiral transition turns into a first order transition, because the 
chiral transition of Np = 3 QCD at m uc i = m s = is expected to be first order. For 
m s < m*, the second order edge of the first order region will deviate from the vertical axis 
according to m u d oc (m* — m s ) 5 / 2 . In phenomenological studies of heavy-ion collisions, 
it is important to know the location of the physical point in Fig. |l|, for which a definite 
prediction is not available from effective theories. 

Identification of the physical point is not easy on the lattice either because it requires 
a precise scale determination at small quark masses in Np = 2 + 1 full QCD. Previous 
studies using staggered quarks suggest that the physical point locates in the crossover 
region [ |ll|] , while a study using Wilson quarks suggests a first order transition [ |12 
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Figure 2. N F = 2 + 1 and 3 staggered quark simulations using improved actions [ [L5|. (a) 
Chiral order parameter at N t = 8 in iVp = 2 + 1 QCD. The bursts are linear extrapolations 
in m U( i to the chiral limit, (b) Triplet quark number susceptibility as a function of N t in 
3-flavor QCD. 



These simulations are done on N t = 4 lattices using the standard lattice actions. Because 
lattice artifacts are expected to be large in these simulations, studies at larger N t using 
improved actions are called for to draw a definite conclusion. 

At the Lattice 2002 conference, several new 2 + 1 flavor simulations were reported. 
The Bielefeld- Swansea Collaboration and the Columbia group presented updates of the 



staggered quark phase diagram at N t = 4 using the standard action [ [13], . The locus 
of the second order line is determined around the Np = 3 critical point using the Binder 
cumulant method (see Fig. [j](b)). The results are consistent with the previous studies 
with standard staggered quarks at N t = 4. 

The MILC Collaboration presented results with staggered quark on N t = 6 and 8 
lattices (12 3 x 6 and 16 3 x 8) using improved actions (1-loop Symanzik gauge action and 
Asqtad staggered quark action) [|I5||. Fixing m s to the physical s quark mass value, three 
light quark masses were studied (m u d/m s = 0.6, 0.4, 0.2). They also performed Np = 3 
simulations at m q /m^ hys = 1, 0.6 and 0.4. Fig. |2](a) shows the quark condensate as a 
function of T and m Bl j. Because of large m u< i, the transition is crossover at all simulation 
points. A crude extrapolation to m u d = shown in Fig. |2|(a) suggests that the transition 
temperature is larger than about 170 MeV there, while, without data at small m U( i, it is 
difficult to conclude about the nature of the transition at the physical point. 

2.2. Toward the continuum limit 

Numerical data obtained on finite lattices must be extrapolated to the continuum limit 
to extract predictions for the real world. In finite temperature QCD, this corresponds to 
the limit of large Nt because the temperature in physical units is given by T = 1/N t a. 
On the other hand, the spatial lattice size N s a should be much larger than T _1 = N t a 
to approximately realize the thermodynamic limit. From experiences in quenched QCD, 
N s /N t ^4 is said to be safe for the EOS. The computation is, therefore, quite expensive 
at large Np This is a reason that major full QCD simulations have been limited to N t 
£6. 
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Figure 3. The EOS from QCQ with 
Wilson-type quarks on N t = 4 (filled sym- 
bols) and 6 (open symbols) lattices using 
improved actions [ M . 
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Figure 4. Continuum extrapolation of the 
energy density at T/T c = 2.5 in SU (3) gauge 
theory on isotropic (£ = 1) and anisotropic 
(f = 2) lattices [ M. 



As mentioned in the previous section, improvements of lattice actions are effective for 
reducing lattice artifacts at small N t . For the EOS, however, improvement seems to be 
insufficient at N t = 4. Fig. |^ shows the energy density in QCD with Wilson-type quarks 
using improved actions (a renormalization-group improved gauge action combined with 
a clover-improved Wilson quark action) [ |4]] . The discrepancies between the results at 
Nt — 4 and 6 are as large as 40-50%, which implies that the N t = 4 values suffer from 
large lattice errors. The discrepancies are too large to attempt a continuum extrapolation 
Nt — > oo using data at Nt = 4 and 6. 

On the other hand, the results obtained on N t = 6 lattices were suggested to be close 
to the continuum limit [ f|, because the EOS from the staggered quark action [ |IBJ are 
quite similar on Nt = 6 lattices. Thus, a precise continuum extrapolation may be possible 
if additional data at Nt = 8 are generated.f] 

The recent study by the MILC Collaboration [0 using improved staggered-type quarks 
seems to be supporting this expectation. Although the EOS has not been calculated, their 
results for triplet quark number susceptibility for Np = 3 clearly show that the N t = 6 
and 8 data are roughly consistent with each other while N t = 4 data clearly deviates from 
them (Fig. 0(b)). 

Calculation of the EOS at Nt > 8 is, however, quite expensive. A possible solution of 
the problem was proposed in Ref . [ 113] . The idea is as follows: a study of thermodynamic 



quantities in the high temperature limit shows that the leading lattice errors from the 
temporal lattice cut-off are larger than those from the spatial cut-off; therefore, the lattice 
artifacts in these quantities may be efficiently reduced on anisotropic lattices with a t < a s . 
From this study, £ = a s /a t = 2 is found to be optimum. 

The idea was tested in quenched QCD [ |THJ. In Fig. [|, the energy densities obtained on 
isotropic (open circles) and anisotropic (filled circles) lattices using the standard plaquette 



2 Here, it is worth to note that a ratio of observables, such as the Wroblewski parameter X s w Xs/Xu 
[ [I?]], can be less sensitive to lattice errors. For these quantities, accurate estimations may be possible 
even at small Nt- 
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action are shown together with that on isotropic lattices using a renormalizat ion-group im- 
proved action (open squares). Results of continuum extrapolations, assuming the leading 
l/iV t 2 scaling law, are also given by corresponding symbols at £/N t = 0. The right-most 
open circle, obtained on an N t = 4 isotropic lattice, shows a large deviation from the 
continuum limit, and, simultaneously, deviates from the leading scaling line. Therefore, 
the continuum extrapolation had to be made using only two data points at N t = 6 and 8 
excluding N t = 4. The filled circles are the results of anisotropic lattices. The advantage 
of using anisotropic lattices is apparent. On the coarsest lattice N t /C, = 4, finite lattice 
spacing error at £ = 2 is much smaller than that at £ = 1 with the same plaquette ac- 
tion. Furthermore, even at Nt/£ = 4, the data is well described by the leading scaling. 
Therefore, we can confidently perform a continuum extrapolation using three data points, 
resulting in smaller and more reliable errors in the continuum limit. Actually, this is the 
first well-controlled continuum extrapolation of the EOS in QCD even in the quenched 
approximation. This anisotropic calculation of the EOS was about factor 5 more efficient 
in computational cost than isotropic calculations. 

Anisotropic lattices may help also in QCD with dynamical quarks. Here, improvement 
of lattice actions is also important. A status report toward this direction was presented 
by Umeda at Lattice 2002 [ 0. 

2.3. Spectral functions by MEM 

Hadronic spectral functions are important in various phenomenological analyses. It is 
desirable to determine them on the lattice directly from the first principles of QCD. In 
Euclidian space-time, the spectral function / is related to a two-point function D by a 
kind of Laplace transformation: 

DM = (0\O(t)O\Q)\0) = / dcuK(uj,T)f(uj), K(u, r) = e~ WT + e~ w{T ~ T \ (1) 

J o 

where O is a hadronic operator and T _1 = N t at- However, because only a finite number 
of discrete data points are available for D, determination of the continuous / from D is 
an ill-posed problem. 

Recently, a new method, the maximum entropy method (MEM), has been started to 
be applied to calculate / from the lattice data. MEM has been successfully applied to 
solve similar ill-posed problems in other fields, such as the image reconstruction problem 
in astrophysics or detective works. 

In MEM, we calculate the "most probable" / by maximizing the probability P[f\DH] 
of / when the data D and prior knowledge H on / are given. By Bayes' theorem, we can 
decompose P[f\DH] as P[f\DH] oc P[D\fH] x P[f\H}. The first part is the probability 
of D given / and H, and can be evaluated by the conventional least \ 2 likelihood function 
L as P[D\fH) oc e~ L [ D 'fl. The second part P[f\H], the probability of / when H is given, 
is something only God knows and concentrates the information what we think "probable" . 
In MEM, the Shannon- Jaynes entropy S is adopted to evaluate 

/•OO 

P[f\H] oc e aS ^ m \ S= du {f(u) - m{u) - f(u) In \f(w)/m(u)]} , (2) 

Jo 

with m(u) > a real function called the default model and a a positive parameter to 
specify the relative weight between L and S. (Note that L attracts f to D while S attracts 
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/ to m.) Eq.(||]) means that we use the term "probable" in the sense of the information 
theory. Here, two additional inputs, a and m(u), are required. The a-dependence turns 
out to be weak in most cases and can be eliminated by an additional statistical argument. 
The m-dependence is also (usually) weak. 

From various tests on the lattice, the method has been shown to work well provided 
that the number of data is large and the statistical quality is high. MEM can correctly 
reproduce the ground state and the first exited state energies from the location of peaks 



in /, as well as the decay constants for these states from the peak area [ |20|, |2T[. When 



D at small distances are included, MEM is accurate enough to reproduce lattice artifacts 
(doublers etc.) too [|21|. It was also demonstrated in a sigma model that MEM is capable 



of describing decay modes [ |2"2 |. A benefit of MEM is that the tasks required to extract 
the information are much lighter than those with the conventional methods. 

It should be nevertheless stressed that, no matter how plausibility arguments are accu- 
mulated, answers to an ill-posed problem are always heuristic. Therefore, it is important 
to understand which properties of / are reliably calculated. For that, let us reexamine 
the role of the input data. From Eq.(|l]), we see that D is essentially the area of / with an 
exponential cut-off at to ~ 1/r. Therefore, we can draw following conclusions: (1) When 
we remove D at t < r min to avoid lattice artifacts, / at to ^c/r min cannot be determined. 



In accord, the choice of m was found to affect / at large to [ [20] . Here, c = 0(1) is a 



constant depending on the precision of D and the required confidence level on /. (A crude 
estimation leads c ~ 2 (8) for the cases of 10% (2%) errors in D at r ~ 1/2T, assuming 
the Born form of / [ |23j.) (2) / at to « 1/T max is also insensitive to the data. Actually, 
when r max is not sufficiently large, we sometimes encounter false peaks at small to [ [H]] . 
See also arguments in [ |23j. (3) The location and the area of a peak will be reliable when 
we input precise data at around r ~ l/to pea ^. (4) To extract more detailed features of 
the peak, such as the width, we need enough number of precise data points within the 
corresponding range of r. 

As stressed in [ [2(J , a rather large number of data points with high statistics is required 
to calculate / precisely. For example, in a study of finite-temperature QCD, N t > 30 was 



argued to be necessary. [ 25|. Because a simulation at large N t is demanding, it is natural 



to adopt anisotropic lattices. In [ |25| , a lattice with £ = 4 was studied. 

Applications of MEM to phenomenological quantities have already been started in the 
quenched approximation of QCD. In [ thermal dilepton rate was calculated from 
a vector spectral function on isotropic Nt = 12 and 16 lattices. From the observed 
suppression of / at small to, they suggested decrease of the dilepton rate relative to the 
Born value at small to (see Fig. |5|). A high number of data points and good statistics seem 
to be critical if a decisive conclusion is to be drawn; however, because a recent higher 
quality data [ ^5] also supports the behavior of / at small to, the decrease of the dilepton 
rate at small to seems to be a real effect. 



3. FINITE DENSITY QCD 

Rich physical content is expected in finite density QCD both theoretically and exper- 
imentally. However, a numerical simulation of QCD at finite chemical potential \x has 
remained to be one of the most challenging issues since the first foundation of lattice 
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Figure 5. Dilepton rate from vector spec- 
tral function in quenched QCD [ |23| . 
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reweighting from T c at // = [ H. 



QCD, because, in the path integral of QCD 

Z = J dUe' Saau9e detM(fx), (3) 

the quark determinant det M becomes complex at /x ^ even in Euclidian space-time: 
det M = | det M\ e ie . 

A trick to validate the statistical treatment of Eq.(^) in Monte Carlo simulations is to 
regard the phase factor e as a part of observables. However, the expectation value of the 
phase factor is found to be exponentially suppressed by the lattice volume, (e w ) oc e~ v . 
Therefore, exponentially large statistics is required as we increase V to keep the accuracy 
of the results. This difficulty is called the phase (sign) problem. 

3.1. QCD at small chemical potentials 

Recently, a (partial) breakthrough of the problem was introduced by Fodor and Katz 
[ ||] applying the reweighting method (spectral density method) to study QCD at small 
densities. Basic idea of the reweighting is simple [ E6[. Because the functional form of the 
measure is known in Eq.fl^), we can correct the distribution of operators measured at a 
point A in the coupling parameter space to that at any other point B as 

dUe S 9 au ge (B) det M , m = dUe -S gauge (A) det M / A x x e S 3auge (A)-S gauge (B) det M ( B \ ( 4 ) 

det M(A) 



target point simulation point 



correction factor 



In practice, however, to obtain a reliable distribution at B, the distributions at A and B 
should have a reasonable overlap. 

Reweighting from \x = (the Glasgow method) is a natural strategy to avoid simulations 
at n 7^ 0. Since early 90 's, however, no attempts made at T = have been able to produce 
a sensible result even at very small \i. We understand that the failure is caused by 
exponentially small width of the distribution in the /i-direction due to the phase problem. 

Here, we can imagine several special cases where we may expect wider distributions: 
(i) Small system volume V. (ii) High temperature, (iii) Phase transition/crossover point. 
These factors will have a cooperative effect. 
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Based on a deep insight of the problem, Fodor and Katz combined all the three and, for 
the first time, succeeded to obtain a reasonable result for QCD at fi ^ [ [2|]. On iV 3 x 4 
lattices with N s = 4, 6 and 8, they performed a simulation at the finite-temperature 
transition point of fi = QCD, and looked for the points with finite overlap in the 
(fi, T) plane as shown in Fig. [| Among these points, they found the endpoint of the 
first order transition line between hadronic and QGP phases at Tg = 160 ± 3.5 MeV and 
He = 725 ± 35 MeV. 

Because the simulation is done with heavy staggered quarks (m^a = 0.025 and m s a = 
0.2) using the standard action, the finite temperature transition is crossover at fi = 0. 
When we decrease the quark masses toward the physical point, the end point of first order 
transition line will move toward smaller \i. As discussed in Sec. [2.1| , unimproved lattice 
actions suffer from sizable lattice artifacts at N t = 4. Whether the first order line reaches 
the T-axis for the physical quark masses is currently unclear. 

In QCD with dynamical quarks, evaluation of the determinants in Eq.([|) is quite time- 
consuming on large lattices. Recently, Allton et al. applied a Taylor expansion method 
to calculate det M at small fi [ |27| , 

bidet MO) = lndetM(O) + Tr (M" 1 ^ ] /!+•••. (5) 




When we adopt the noise method to calculate the quark traces, large lattice volume is 
not a serious problem at least in this part of the calculation. Accordingly, they performed 
a simulation on a 16 3 x 4 lattice in 2 and 3-flavor QCD with improved staggered quarks. 
They confirmed that the method applies well to the RHIC region (fi/T m 0.1), and 
calculated several thermodynamic quantities. In particular, the effect of /i in the EOS is 
estimated to be 0(1%) at the RHIC point. 

Another approach is to try an analytic continuation from pure imaginary fi where det M 
is real [ |28| (see also [ |29|, |30|). This method is also limited to small real \i because (i) 
|Im/i| < 7rT/3 due to a Z(3) periodicity in the imaginary fi direction and (ii) the number 
and statistical accuracy of data are limited. Concretely, we first fit the data at several 
values of imaginary fi by a low-order polynomial in fi, and replace fi by —ifi in the final 
fit-function. Results on 6 3 x 4 and 8 3 x 4 for Np = 2 QCD are reported to be consistent 
with those from other groups. 

3.2. Large chemical potentials 



QCD was suggested to have rich phase structures at high densities [ pl| . The region in 
the phase diagram may be studied by GSI and JHF. The methods described in Sec. [3T| 
however, cannot be applied to study the high-density regions. To obtain insights for QCD 
at large fi, a special case of QCD and several QCD-like theories are studied in which det M 
is real even at fi ^ 0. 

When we consider QCD with isovector chemical potential, fi u = —fid, the functional 
measure [det M(fi) det M(— fi)] is real because det M(—fi) = (det M(/i))*. Simulations in- 
dicate pieces of evidence for the expected phase transition to a phase with 7r + -condensation. 
Results of observables are compatible with predictions from an effective chiral theory [ 
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The quark determinant is also real at fi ^ in two-color QCD [ [33|, gj, g5j, p6| . From 
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simulations, existence of di-quark condensation has been demonstrated and consistency 
with effective chiral theory has been reported. Muroya et al. studied meson masses at 
7^ and showed that m p decreases with increasing \i [ |36|, which may explain the 



dilepton enhancement observed by CERES [ [!7|] if the behavior remains in real QCD. 

Finally, a simulation of the NJL model showed an evidence of a BCS transition in accord 
with the color-superconducting scenario [ BS1. See also [ B^] for other study of NJL. 



4. CONCLUSIONS 

Since the QM 2001 conference, a couple of steps have been made in lattice QCD to 
reduce the distance between theory and experiment. For several quantities, lattice studies 
of QCD at finite temperatures are entering the stage of precise quantitative predictions. 

New 2+1 flavor simulations have further constrained the finite-temperature QCD tran- 
sition, although the clarification of the nature of the transition at the physical point 
has remained an open problem. Simulations at large N t suggest that a precise contin- 
uum extrapolation of thermodynamic quantities may be possible with data at Nt = 6, 8, 
■ ■ •. Calculation of EOS at N t > 8 is, however, demanding. A way out will be given by 
anisotropic lattices, which have enabled us to carry out the first well-controlled continuum 
extrapolation of the EOS in quenched QCD. An extension to QCD with dynamical quarks 
is currently being attempted. The MEM will enlarge the predictive power of lattice when 
used carefully, and the first applications to the quark matter phenomenology have already 
begun. The hard numbers, such as T c (7Vf=2) « 170-175 MeV, e {Np=2) /T 4 | T/Tc=1 . 5 « 14, 
e^ Np=2 ^ /T 4 \t/t c =i.5 ~ 3, are unchanged through these developments, but the confidence 
in these values has become much firmer. 

In the last decades, the finite density QCD has been a barren area in lattice QCD in spite 
of its great importance and big efforts, due to a technical difficulty of numerical simulations 
at finite \x. Last two years, however, a (partial) breakthrough has been made for the case 
of small chemical potentials around the finite temperature transition point, and the first 
sensible results for QCD at ji ^ are obtained. Different groups, adopting different 
methods, have reached a general agreement of results. Further studies are underway to 
reduce the lattice artifacts. These methods cannot be simply applied for large //. To 
obtain insights for QCD at large /i, several QCD-like theories, for which the conventional 
technique is applicable also at large //, are being studied. Results obtained so far seem to 
be consistent with the color-superconducting scenario. 
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valuable discussions, comments and data. This work is supported in part by Grants-in-Aid 
of the Ministry of Education (Nos. 12304011 and 13640260). 
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